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ABSTRACT 



m 

o 

^D I In the continuum the Bianchi identity imphes a relationship between different components of the 

curvature tensor, thus ensuring the internal consistency of the gravitational field equations. In this 
CT^i paper the exact form for the Bianchi identity in Regge's discrete formulation of gravity is derived, 

KTv' ^y considering appropriate products of rotation matrices constructed around null-homotopic paths. 

The discrete Bianchi identity implies an algebraic relationship between deficit angles belonging 

^ ■ to neighboring hinges. As in the continuum, the derived identity is valid for arbitrarily curved 

H ■ 

C^ I manifolds without a restriction to the weak field small curvature limit, but is in general not linear 

in the curvatures. 



(Submitted to Class, and Quant. Gravity) 



^e-mail address : hhamber@uci.edu 
e-mail address : gkagel@uci.edu 



1 Introduction 

In this paper we investigate the form of the Bianchi identities in Regge's [1] lattice formulation 
of gravity [2, 3, 4, 5, 6, 7, 8, 9, 10]. The Bianchi identities play an important role in the contin- 
uum formulation of gravity, both classical and quantum-mechanical, giving rise to a differential 
relationship between different components of the curvature tensor. It is known that these simply 
follow from the definition of the Riemann tensor in terms of the afhne connection and the metric 
components, and help ensure the consistency of the gravitational field equations in the presence 
of matter. At the same time they can be regarded as a direct consequence of the local gauge 
(diffeomorphism) invariance of the gravitational action, since they can be derived by invoking the 
invariance of the action under infinitesimal gauge transformations (see for example [18, 19, 20, 21]). 

In this paper we will show that the lattice formulation of gravity has an equivalent form of 
the Bianchi identities, which are both exact, in the sense that they are valid for arbitrarily curved 
lattices, and reduce to their continuum counterparts in the weak field limit. We will derive the 
exact lattice Bianchi identities in three (Section 4) and four dimensions (Section 5) explicitly 
by considering the product of rotation matrices along paths which are topologically trivial (i.e. 
reducible to a point). By expressing the rotation matrices about each hinge in terms of the local 
representatives of the curvatures, namely the deficit angles, we will obtain an algebraic relationship 
between deficit angles, area and volumes pertaining to neighboring simplices. 

For lattices which are close to fiat, it will be shown that the derived set of identities is analogous 
to the Bianchi identities in the continuum, once the edge lengths are identified with appropriate 
components of the metric in the continuum (Sections 9 and 10). We will therefore extend and 
complete previous results on the lattice Bianchi identities, which so far have been restricted to the 
weak field limit [1, 3, 4, 22, 23]. The results presented in this paper should therefore be relevant for 
both classical (see for example [24, 25, 26, 27] an references therein) and quantum (see for example 
[28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38] and references therein) discrete gravity. ^ 

It is well known that in the continuum the Bianchi identity for the curvature tensor ensures the 
consistency of the Einstein field equations. For the Riemann curvature tensor the un-contracted 
Bianchi identities read 

^ ua/3;f ~^ -^ U'ya;f3~^ -^ u/3r,a = 0, (1.1) 

^Recent reviews of Regge gravity can be found in [10, 12, 13, 14, 15, 16], while a comprehensive collection of up-to- 
date references is assembled in [17]. Some further mathematical aspects of piecewise linear spaces, with some relevance 
to lattice gravity, are discussed in the above cited references [5, 6], as well as in [39, 40, 41, 42, 43, 44, 45, 46, 47] and 
references therein, and more recently in [48, 49, 50, 51, 52, 53] and references therein. 



or more concisely 



-^ iy[al3; 7] ~ ^ ' 



(1.2) 



where [• • •] denotes symmetrization. These identities are easily derived by inserting into the above 
expression the explicit definition for the curvature tensor in terms of the metric ^f^jy, 

1 d^gx. 



It is then easy to see that in d dimensions there are 

d{d-l) /d\ _ d'^{d-l){d-2) 
2 [3) ~ 12 

un-contracted Bianchi identities, and thus 3 identities in d = 3 and 24 identities in d = 4. In their 
contracted form, the Bianchi identities imply for the Ricci tensor 



(1.3) 



(1.4) 



J^va; 7 -tiuy; a ~T~ -ti ^^q; ^ 



0, 



and for the scalar curvature 



/l;^ 2K ^. ^ 



0. 



These relations in turn give the contracted Bianchi identity 



i2^-in^ 



0, 



(1.5) 



(1.6) 



(1.7) 



which always corresponds to d equations in d dimensions. We note here that a simple physical 
interpretation for the Bianchi identity can be given in terms of a divergence of suitably defined 
stresses [18]. Thus for example in three dimensions one has 



pM 



with P^ 






[l.i 



It is also well known that the Bianchi identities are required for ensuring the consistency of the 
gravitational field equations. Consider the classical field equations with a cosmological constant 
term, 



R 



pu 



\gpvR + J^g^u = SirGT^^, 



(U 



with A = 8ttGX the cosmological constant. Applying a covariant derivative on both sides one has 



R'^.-^S'^^R 



S^GT^, 



v; fii 



(1.10) 



which for a covariantly conserved energy-momentum tensor 



T^ =0 



(1.11) 



is only consistent if the contracted Bianchi identity of Eq. (1.7) is identically satisfied. 

In d dimensions one has d contracted Bianchi identities. Since there are in general d{d + l)/2 
equations of motion, as well as d harmonic gauge fixing conditions, one has for the number of 
independent gravitational degrees of freedom in d dimensions 

*l±i) -„-<,= *L:1 (1.12) 

2 2 ^ ' 

which indeed reproduces correctly in four dimensions the two independent helicity states appropri- 
ate for a massless spin two particle. 

There exists also a close relationship between the Bianchi identity and the gauge invariance 
of the gravitational action. We shall take note here of the fact that the Bianchi identity can be 
derived from the requirement that the gravitational action 

1 



Ig[9] = -j^;^] d^xV^Rix) , (1.13) 

being a scalar, should be invariant under infinitesimal local gauge transformations at a space-time 
point X, 

^9,iu{x) = -g^i\{x) dyX^{x) - g\u{x) d^x^{x) - dxg^y{x) x^{x) . (1.14) 

After substituting the above expression for a gauge deformation into the variation of the action 
given in Eq. (1.13), 

5Ig[9] = ji^Jd'xVaix) [R'^'^-lg^^^R] 6g,, , (1.15) 

one obtains again, after integrating by parts, the contracted Bianchi identity of Eq. (1.7). 

Let us now turn to the lattice theory. In a d-dimensional piecewise linear space-time the 
expression analogous to the Einstein-Hilbert action was given by Regge [1] as 

Ir= Y. ^t'^ ^h , (1.16) 

hinges h 

where Af^ is the volume of the hinge and dh is the deficit angle there. The above lattice action 
is supposed to be equivalent to the continuum expression 

Ie = \j d'^x^R , (1.17) 

and indeed it has been shown [3, 5, 11] that Ir tends to the continuum expression as the simplicial 
block size (or some suitable average edge length) tends to zero in the appropriate way. The above 



Regge form for the lattice action can be naturally extended to include cosmological and curvature 
squared terms [7, 10] 

1(1^) = J2[^^h-k^h6h + a6lAl/Vh + --] , (1.18) 

h 
with k^^ = 8ttG. In the limit of small fluctuations around a smooth background, /(/^) corresponds 
to the continuum action 

I[g] = Jd^x^[X-^R + ^ R^upaR^"'" + ■■] ■ (1.19) 

In the following we will focus on the Regge term (proportional to k) only. 

Variations of Ir in Eq. (1.16) with respect to the edge lengths then give the simplicial analogs 
of Einstein's equations, whose derivation is significantly simplified by the fact that the variation of 
the deficit angle is known to be zero in any dimensions, 

6Ir = Y.^{At'^)5n , (1.20) 

h 

as happens in the continuum as well (where one finds that the variation of the curvature reduces 
to a total derivative). In three dimensions the above action gives for the equation of motion 6h = 
for every hinge in the lattice, whereas in four dimensions variation with respect to Ip yields [1] 

^lpJ2 5h cot 9ph = (1.21) 

hDlp 

where the sum is over hinges (triangles in four dimensions) labeled by h meeting on the common 
edge p, and 9ph is the angle in the hinge h opposite to the edge p. The above equation is the lattice 
analog of the field equations of Eq. (1.9), for pure gravity and vanishing cosmological constant. 

Numerical solutions to the lattice equations of motion can in general be found by appropriately 
adjusting the edge lengths according to Eq. (1.21). Since the resulting equations are non-linear in 
the edge lengths, slight complications can arise such as the existence of multiple solutions, although 
for sufficiently weak fields one would expect the same level of degeneracies as in the continuum [3] . 
Several authors have discussed the application of the Regge equations to strong field problems in 
classical general relativity, and some samples can be found in [24, 25, 26]. The relevance of the 
Bianchi identities to a numerical solution to the lattice field equations - using for example a 3 + 1 
time evolution scheme - resides in the fact that they are in principle a powerful tool to check the 
overall accuracy and consistency of the numerical solutions. 

The Bianchi identities also play an important role in the quantum formulation. In a quantum- 
mechanical theory of gravity the starting point is a suitable definition of the discrete Feynman path 



integral [8, 7, 10]. In the simplicial lattice approach one starts from the discretized Euclidean Feyn- 
man path integral for pure gravity, with the squared edge lengths taken as fundamental variables, 

^L = r U (.Vd{s)r n '^^^J 0['|] exp I - ^ (a K - A: 6hAh + a SIaI/Vh + •••)[ . (1.22) 

The above regularized lattice expression should be compared to the continuum Euclidean path 
integral for pure gravity 

Zc = III {^gi^y n d9f.u{x) e^p{- fd'x^[X-^R+JR^.,.R^^'''^'' + ■ 
(1.23) 
In the discrete case the integration over metrics is replaced by integrals over the elementary lattice 
degrees of freedom, the squared edge lengths. The discrete gravitational measure in Zl can be 
considered as the lattice analog of the DeWitt continuum functional measure [5, 6, 8, 9, 12, 38]. 
A cosmological constant term is needed for convergence of the path integral, while the curvature 
squared term allows one to control the fluctuations in the curvature [7, 10]. In the end one is 
mostly interested in the limit a ^ 0, where the theory, in the absence of matter and after a suitable 
rescaling of the metric, only depends on one bare parameter, the dimensionless coupling fc^/A. 

In the quantum theory the Bianchi identities of Eq. (1.7) are still satisfied as operator equations, 
and ensure the consistency of the quantum equations of motions. In ordinary lattice nonabelian 
gauge theories an attempt has been made to entirely replace the functional integration over the 
gauge fields by an integration over field strengths, but now subject to the Bianchi identity constraint 
[54]. In the case of gravitation such an approach is more difficult, since the analog of the gauge 
field is represented by the affine connection, and not the by the curvature tensor. 

2 Lattice Parallel Transport 

To construct the lattice Bianchi identities we will follow a strategy similar to the one used in 
the derivation of the exact lattice Bianchi identities in non-abelian lattice gauge theories. There 
the Bianchi identities can be obtained by considering the path-ordered product of SU{n) gauge 
group rotation matrices, taken along a suitable closed path encircling a cube. The path has to be 
chosen topologically trivial, in the sense that it can be shrunk to a point without entangling any 
plaquettes [54, 55]. 

Let us therefore first review the notion of parallel transport of a test vector around a small loop 
embedded in the lattice. Consider a closed path T encircling a hinge h and passing through each 
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of the simplices that meet at that hinge. In particular one may take T to be the boundary of the 
polyhedral dual area surrounding the hinge. For each neighboring pair of simplices j,j + 1, one 
can write down a Lorentz transformation L^, which describes how a given vector (p^ transforms 
between the local coordinate systems in these two simplices, 

<P'^ = [L(j,i + l)j V • (2.1) 

Now in general it is possible to choose coordinates so that L^ is the identity matrix for one pair 
of simplices, but then it will not be unity for other pairs. The above Lorentz transformation 
is directly related to the continuum path-ordered (P) exponential of the integral of the affine 



connection (r> 



fiX 



by 



L. 



J path r^O 

Pp between simplices 



(2.2) 



The connection here is intended to only have support on the common interface between the two 
simplices. 

Next we will consider the product of rotation matrices along a closed loop T. The path can 
entangle several hinges, or just one, in which case it will be called a closed elementary loop. On 
the lattice the effect of parallel transport around a closed elementary loop T is obtained from the 
matrix [7] 



[l[LiJ,J + l] 



flU 



JhUW 



flU 



(2.3) 



rW ;. 



where Uij,u is a bivector orthogonal to the hinge h, defined in four dimensions by 

1 



U^H) 



irj7.. 



^^ 2 Ah ^"P"" {o.y{b) ' 



(2.4) 



with /f N and Z?^-, two vectors forming two sides of the hinge h. We note that in general the validity 
of the lattice parallel transport formula given above is not restricted to small deficit angles. For a 
closed path F, the total change in a vector i/)^ which undergoes parallel transport around the path 
is given by 



+ 



n Mi,i + i: 



(2.5) 



pairs of 
simplices on F 



For smooth enough manifolds, the product of Lorentz transformations around a closed elementary 
loop F can be deduced from the components of the Riemann tensor, 

(2.6) 



n L{j,j+i) 


U 


' (.Rpa^'"' 


pairs or 




simplices on T 







where (R' p^Ya — ^^npa i^ ^^^ curvature tensor and S'"^ is a bivector in the plane of F, with 
magnitude equal to 1/2 times the area of the loop T. (For a parallelogram with edges a'' and h^^ 
S^"" = 2(0'^^'' — a^b^)). The above result then reproduces to lowest order the parallel transport 
formula 

5<p^ = i?;,.S'"^0. . (2.7) 

Comparison of Eq. (2.3) and Eq. (2.6) means that for one hinge one may make the identification 

R^.upa^"' ^ ShUJ!:;^ . (2.8) 

It is important to notice that this relation does not give complete information about the Riemann 
tensor, but only about its projection in the plane of the loop T, orthogonal to the given hinge. 
Indeed the deficit angle divided by the area of the loop can be taken as a definition of the local 
sectional curvature K^ [10] 

"h _ T^ _ n,pupcr t^gCfeCgCfc /2 q\ 

^Fh (OppQua ~ Qfiagupj^ae^eae^ 

which represents the projection of the Riemann curvature in the direction of the bivector e^ A eb. 
The lattice Bianchi identities are derived by considering closed paths that can be shrunk to a 
point without entangling any hinge. The product of rotation matrices associated with the path 
then has to give the identity matrix [1, 3]. Thus, for example, the ordered product of rotation 
matrices associated with the triangles meeting on a given edge has to give one, since a path can be 
constructed which sequentially encircles all the triangles and is topologically trivial 

= 1 . (2.10) 

Other identities might be derived by considering paths that encircle hinges meeting on one point. 

3 Geometric Setup 

The discrete analogue of the Bianchi identity will be derived by considering a product of ro- 
tation matrices along a homotypically trivial path. This section discusses the general geometric 
setups needed to define correctly the product of rotation matrices entering the exact lattice Bianchi 
identities derived later in the paper (Sections 4 and 5). For the 3-d case, consider a tetrahedron 
with a point in its interior. Connect the vertices of the tetrahedron to the point in the center. We 



n 


'e^nui^y 


hinges h 


meeting on edge p 





now have formed 4 tetrahedra from the original tetrahedron. In three dimensions, hinges are edges, 
so here we have enclosed four hinges: one connecting each vertex of the original tetrahedra to the 
interi' 




Fig. 1. In three dimensions four tetrahedra meet on a point, labeled here by 0. Deficit angles 5{)i, 5^2, <^03 
and (5o4 are associated with edges 0—1, — 2, — 3 and — 4 respectively. The Bianchi identities are obtained 
by taking an ordered product of rotation matrices along a path which encircles all four hinges (edges here) 
and is topologically trivial, in the sense that it can be shrunk to a point. 



Referring to Fig. 1, for the moment only considering flat space, call the interior point (zero) 
and place it at the origin. Let us take our coordinate system so that vertex 1 lies on the positive 
z axis (zl > 0). Let us take vertex 2 to lie in the x-z plane with z2 < and x2 > 0. Let us take 
vertex 3 to have zS < 0, xS < and y3 > 0. Finally, let us take vertex 4 to have z4 < 0, x4 < 
and 2/4 < 0. So, to summarize: 



Vertex : 1 2 3 4 

^'' n ~ ~ ~ (3-1) 

X : + - - ^ ' 

y: 0+- 

Now, the following argument will apply for all cases where we have the center point completely 

surrounded; the aforementioned restrictions in Eq. (3.1) are mentioned only to give the reader a 

nice picture of the situation. 

Now, "curve the space". This is done by changing one of the edge lengths. Any 9 of the 10 
edge lengths can be chosen arbitrarily (provided the center point is completely surrounded by the 
constructed volumes and provided that real areas and real volumes are still formed) and the space 
will still be flat. So curvature, in this setup, just amounts to adjusting one edge length {I34, the 
edge between vertices 3 and 4, is the easiest one to adjust). For an arbitrary setup, one adds 
more edges until the relative flat space locations of all vertices are specified, making sure not to 
add any edges between points whose relative flat space location is already determined but rather 
marking each such edge as an "unadded edge" (and making sure that all areas and volumes are 
real); then, all remaining ("unadded") edges are determined for flat space, and it is the varying of 
those remaining edges which is the source of all curvature for that setup. This last comment applies 
to d-dimensions and general lattices, not just d-simplices surrounding one point, and addresses the 
question of curvature invariance under edge variation, i.e. curvature invariance under the edge 
variation of the added edges (with the "unadded" edges still unadded) . 

Further, let us label the four tetrahedra by a point in each of them along the line const x (li + 
I2 + I3) where li, I2 and I3 are vectors based at the point 0, the interior point. 
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'BDC 



Fig. 2. A topologically trivial path which encircles four hinges in three dimensions, and can be shrunk to a 
point. The vertices A, B, C and D reside in the dual lattice, and can conveniently be placed at the centers of 
the tetrahedra shown in Fig. 1. It is then advantageous to label the deficit angles by the vertices in the dual 
lattice. The path ordered product of rotation matrices around the shown path then reduces to the identity 
matrix. 



Referring to Fig. 1, and Fig. 2 which, with triangle BCD behind A, views Fig. 1 from 
underneath, let us call the point in the center of the tetrahedra formed via vertices 0123, and not 
4, point D; similarly call the point at the center of 0234, and not 1, point A, that at the center 
of 0134, and not 2, point B and that at the center of 0124, and not 3, point C. We also have 
tetrahedron A which is the tetrahedron containing point A, tetrahedron B containing point B, etc. 
These four points, point A, point B, point C and point D form a Voronoi tetrahedra. Going from 
A ^ B ^ C ^ A, for example, goes around hinge 04. In general 

hinge gone around rotation Voronoi vertex not in path positive path negative path 

01 m A B ^C ^ D^ B B ^ D ^C ^ B 

C^ D^A 
D^ B ^ A 
B^C ^ A 
(3.2) 

so that rotations which are clockwise when viewed from outside our setup, or counterclockwise 

when viewed from the point in the middle, are associated with a "-R", and rotations which are 

counter-clockwise when viewed from outside our setup, or clockwise when viewed from the point 

in the middle, are associated with a "+R", as is traditional for a right handed coordinate system 

inside the set-up. 



02 


R2 


03 


m 


04 


R4 
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A- 


^ D- 


-^C - 


^ A 


A 


C 


A- 


^ B- 


^ D- 


^ A 


A 


D 


A- 


^C - 


^B - 


^ A 


A 
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ADB 



'ABC 



B 



BDC 



Fig. 3. Paths associated with different hinges in three dimensions are shown in different shades of gray. The 
product of rotation matrices around a single hinge is simply elated to the deficit angle for that hinge. In 
order to obtain a closed path, segment AB has to he transversed more than two times. 



Now, let us try and compose a null path around all four hinges. Referring to Fig. 3, we write 
our null path as 

A^B^D^A,A^D^C^A,A^C^B^A,A^B,B^C^D^B,B^A (3.3) 

or, noting that the path is a null path, we have 

1 = Rtot = T^Ib R^b Ta-^b R4: R2 R3 (3.4) 

where Ta^b is a rotation matrix representing the rotation that occurs when a vector is parallel 
transported from tetrahedron A directly to tetrahedron B, Rtot means the total rotation matrix 
after going through the whole path, and RIb is Rl written in B's coordinate system; the other 
rotations, with no subscripts, are written in A's coordinate system. Eq. (3.3) can be written 
completely in A's coordinate system as 



1 = R^^^ = Rl Ri R2 R3 



(3.5) 



Also in reference to the F's, there is a subtlety to note. The F's are always "direct" F's between 
adjacent d-simplices. For example, F^^^ is not Ta^c^b (the latter being starting at A, then go 
directly to C, and then go directly to B); indeed, since 



Rl = Ta-^bT-/ 



A^C^B 



(3.6) 
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both cannot simultaneously be set equal to the unit matrix. 

In general, when setting up a d-dimensional global coordinate system, one specifies the coordi- 
nate system in one d-simplex and N^ — 1 F's (Nd is the number of d-simplices in the lattice) such 
that all d-simplices are linked directly or indirectly; to specify any more F's would be to specify a 
pre-determined (via edge lengths) loop rotation. One can, for e.g., choose all N^ — 1 F's equal to 
the unit matrix, so that Tj^^b in Eq. (3.4) would be one of three F's set equal to the unit matrix. 

If one only specifies F's between adjacent d-simplices, one can think of this coordinate system 
in fiat d-space as d-simplices which 1) are attached to each other only by the (d-l)-simplex "faces" 
corresponding to the N^ — 1 F's, and 2) can go through each other (for e.g., consider a hinge with 
greater than 2tt radians around it) . Of course, when viewed in the actual curved space, the d- 
simplices do not go though each other. Some sample coordinate systems would be obtained by 1) 
specifying F's between one specific tetrahedron (or d-simplex) and every other one, such as Tj^^b, 
F^^C* and Ta^d in our particular set up, and 2) specifying the F's in a chainlike fashion, such as 
Ta^b, r_B-»C) and Tq^d in our set up. 




13 



Fig. 4- Similar to Fig. 1., but now in four dimensions. In Jour dimensions five simplices meet on a point, 
labeled here by Q . Deficit angles 5qi2 etc. are associated with triangles 0—1 — 2 etc. respectively. The Bianchi 
identities are obtained by taking an ordered product of rotation matrices along a path which encircles several 
hinges (triangles here) and is topologically trivial, in the sense that it can be shrunk to a point. 

Now, let us consider one four-dimensional set up in particular. Consider a four simplex with a 
point in the middle. This divides the original 4-simplex into 5 new 4-simplices. Then, readjust the 
edge lengths and curve the space (see Fig. 4). We now consider parallel transporting the vector 
around a null path within this simplicial complex. Label each of the new 4-simplices A, B, C, D 
and E (see Fig. 5). 




Fig. 5. Similar to Fig. 2., but now in four dimensions. The topologically trivial path which encircles all 
shown hinges can be shrunk to a point. The vertices A, B , C , D and E reside in the dual lattice, and can 
conveniently be placed at the centers of the simplices shown in Fig. 4-- As in the three-dimensional case, it 
is advantageous to label the deficit angles by the vertices in the dual lattice. The ordered product of rotation 
matrices around the shown path then reduces to the identity matrix. 

We now make an important point for this particular set up: the number of 4-simplices around 
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a triangle hinge is 3. There are 3 vertices in the two dimensional space of the hinge and three 
vertices in the remaining two dimensions. The three vertices in the hinge and any two of the three 
vertices outside the hinge form a 4-simplex, so that three 4-simplices surround the hinge. ^ 




ABD 



'BCD 



Fig. 6. Similar to Fig. 3., hut now in Jour dimensions. Paths associated with different hinges are shown 
here in different shades of gray. The product of rotation matrices around a single hinge is simply related to 
the deficit angle for that hinge. Again, in order to obtain a closed path segment AB has to he transversed 
more than two times. In constructing a product of rotations, it is sufficient to consider only four of the five 
vertices; here, vertex E is excluded. In higher dimensional analogues, one has vertices F, G, etc., which are 
also not used. 

So, we now consider a path which includes 4 of the 5 simplices which surround the "point in 
the middle". Specifically, referring to Fig. 6, the path is 

A^B^D^A,A^D^C^A,A^C^B^A,A^B,B^C^D^B,B^A (3.7) 

The same path works for 3-d, which then includes all 4 of the 4 tetrahedra surrounding the point 
in the middle (see Eq. (3.3)). 



^This three d-simplex result is now easily seen to be true in any number of dimensions for the point in the middle 
set up. 



15 



4 Product of Rotations in Three Dimensions 



Having described the general geometric setup in the previous section, we now proceed to give 
an exphcit form for the rotation matrices in three dimensions, and derive an exact form of the 
lattice Bianchi identity by considering a product of rotation matrices along a homotypically trivial 
path. It turns out to be very convenient to be able to express a rotated vector in terms of the "old" 
vector and the hinge edge. One notes that the only rotation occurs to the portion of the vector 
that is perpendicular to the hinge, so that (v- 1)1 is part of the new vector. One proceeds to form an 
orthogonal coordinate system with the old vector, -u, and the hinge, /, using l,v x / and I x {v x I). 
So, since v clearly has no component parallel to v x I, the component of v which is rotated in the 



V X I 



I X (v X I) plane is parallel to I x (v x I). So, one finds 



if = (v-l)l + 
which simplifies to 



V ■ [i X {v X i)] 
\ix{vxl)\ 



cos 6 



I X {v X I) 
_\ix ivxi)\_ 



sin 5 



V X I 

\v X i\ 



v' 



.2^. 



2sm^-(v-l) 



(4.1) 



(4.2) 



I + {cos6)v + {sm6)l x v 

The total rotated vector (after all successive rotations) is most easily found using a recursive 
application of Eq. (4.2) via standard dot/cross product rules, as well as being sure to set up one's 
coordinate system where the Ty^^s in Eq. (3.4) equals 1. Then, one notes that the most general 
form of a rotation can be written as 



Rtotv = av + b X v + ^(cj • v)di 



a + Ciidii -63 + Ci2dii 62 + Cizdii 

h + Ciidi2 a + Ci2di2 -61 + Ci3di2 I v 

-62 + Ciidis 61 + Ci2di3 a + Qsdis 



(4.3) 



where a, b, and the q are easily found with the expression one had gotten for v'r,-^^^ using Eq. (4.2). 
Here, q runs over all possible vectors in the rotation, which for the sample case we are considering, 
is just li, 12,13, I4 and their 6 cross products. The di are then to be viewed as the vector coefficients 
of the 4 + 6 = 10 Cj-tT's in Eq. (4.3) (look at the non-matrix expression). The a, the b and the di are 
extremely complicated expressions in terms of deficit angles, internal angles, normalized volumes 
and normalized areas. Here, normalization means dividing by the product of the magnitudes of 
the defining vectors. Do note that due to the vector nature of the cross product in 3-d, a does not 
contain normalized areas. 

From the matrix form of Eq. (4.3) one can find a scalar and vector identity. We find these 
equations by noting 

Rtot -1 = (4.4) 
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Taking the trace of both sides while using Eq. (4.3) gives the scalar equation 

3(a-l) + Y,Ci-di = (4.5) 

i 

while taking the antisymmetric part of Eq. (4.3) gives the vector equation 

-fe-^E^^xrf'^^ = (4.6) 

i 

and taking the symmetric part of Eq. (4.3) gives a tensor equation which is unnecessary to discuss. 
So, here, we have found a vector and scalar analogy of the product of rotations for arbitrary 
deficit angles. These (un-contracted) identities are in terms of deficit angles, internal triangle 
angles, normalized areas and normalized volumes. The normalized areas, naturally vectors from a 
geometric point of view, appear explicitly only in the vector identity. Again, it should be noted that 
Eq. (4.4) is only valid for a product of rotations around hinges which gives no rotation. Permute 
two of the rotation matrices making up the total rotation, and the new total rotation matrix will 
not be the unit matrix. 

The explicit form of Rtot in Eq. (4.3) can then be used, in combination with Eq. (4.6), to write 
the 3-dimensional completely contracted Bianchi identities as 



which simplifies to 



J^^{2b + J2c,xdi)^ = (4.8) 



Let us now consider the small deficit angle limit to first order. This simplifies Eq. (4.2) dramatically 
to 

Kmi,\\5 = V + ^{i X v) (4.9) 

which gives, after applying successive rotations, 

a = 1 

b = EiSJi (4.10) 

c = 

so that Eq. (4.6) becomes 

-26 = (4.11) 

which, following a similar derivation to that of Eq. (4.8), allows the fully contracted Bianchi iden- 
tities to be written as 

^ b^ = (4.12) 



Viv) 
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Eq. (4.5) is a trivial 0=0 to first order. To second order, it becomes 

-b-b = (4.13) 

where the b used is just the h in Eq. (4.10). This equation can clearly be seen as a consequence of 
Eq. (4.11), a first order equation. Eqs. (4.11) and (4.12) (as well as Eq. (8.10), reduced to three 
dimenions) have been verified to 1st order, and Eq. (4.13) has been verified to 2nd order. But for 
higher orders, these equations have been shown to be violated in favor of, repectively, Eqs. (4.6), 
(4.8), and (4.5) (as well as Eq. (8.11), again reduced to three dimensions) which are valid in all 
cases. Additionally, all these equations are true for arbitrary set-ups, not just our point in the 
middle set-up. 

5 Product of Rotations in Four Dimensions 

The four-dimensional case is quite similar to the three-dimensional case discussed in the previous 
section, and no major additional complexities arise. We will give here again an explicit form for 
the rotation matrices in four dimensions, and derive an exact form of the four-dimensional lattice 
Bianchi identity by considering a product of rotation matrices along a homotypically trivial path. In 
the four dimensional case, considering how a vector rotates when it is parallel transported around a 
hinge, we can, as in the three dimensional case, form again an orthogonal coordinate system using 
the (old) vector and the hinge, the additional edge in the hinge compensating for going up one 
dimension. Taking li and I2 as the two edges which form the hinge, one uses the following four 
vectors to form the orthogonal coordinate system: li, I2, v y.lixl'2 and /2 x ^1 x ("W x '1 x ^2) where 

,7 h — {h • h)h h — {h-h)h /^ ,n 

io = — — - ^ = — (5.1) 

\h-{h-h)h\ 2A[hM] 

and where A[/i, l2\ is the area of the triangle, not the parallelogram, formed by /i and I2. 

Now, in four dimensions, the hinge we rotate about is a triangle, and the plane in which the 
rotation occurs is perpendicular to that triangle. Since li and ^2 span the space of the triangle, 
the components of the rotated vector in those two perpendicular directions will not change. Now, 
since the only remaining component of the vector is in the ^2 x /i x (iJ x /i x /2) direction, v clearly 
being perpendicular to v xli x I2, we can write 
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,^c-sf /^r/Nf/ f ^- [fo X /i X (w X fi X /o)] 1 f 
(7;-/i)/i+(t;4)/^+<^ —^ -^ ^ \lcos5 

y |i2 X ^1 X (w X /i X igjl J t 



which simphfies to 



/2 X /i X (tJ X /i X ^2) 
1/2 X /i X (fJ X /i X /2)| 



+ sin (5 



V xlixl'2 

\v xlix /2I 
(5.2) 



tJ' = 2 sin^ -[(■«• /i)fi + (w • /I)/!] + [cos(5] ^+ [sin (5] -u x / 1 x i'2 (5.3) 

Now, as in the three dimensional case, we consider null paths as being products of rotations as we 
parallel transport a vector; only now we are going through 4-simplices as opposed to tetrahedra. 
Recall that in the three dimensional case, the path could be written in terms of the four tetrahedra 
it went through, and only depended on the fact that each tetrahedron was directly connected to 
every other tetrahedron. We can follow the same procedure in the four dimensional case by choosing 
four of the five 4-simplices and applying Eq. (3.7) (we could interchange 4-simplex E with any of 
the other 4-simplices in Eq. (3.7)). We have 

4 
Rh,tot -1= Y{RH,m-I = (5.4) 

m=l 

where the li is the same as in the previous equations. In general, the li index can be taken to 
represent the "hinge base" for the product of rotations, i.e. the edge contained in all hinges involved 
in the product of rotations; the 1'2 of the previous equations, on the other hand, varies from hinge 
to hinge. In our set-up, li is the edge "opposite" the excluded 4-simplex. (To be sure that one does 
not have the order of rotations for rotating by +(5's reversed, one can either project out along li or 
simply try both orders). 

Following the 3-d case, we write the most general product of matrices, and hence of rotation 
matrices, as 



[jklmax 

n i k ifc=l i 



( a + c\d] Hk.ibikMjf^ + cUl Efc,(&ifc,&2jr + cM Y.kj{hkM3?^ + 4dl\ 

£(5.5) 



Efc,(&ifc, h2,f^ + c]d1 a + 4d} Ek.ihk, b2jy' + cf df Ekjihk, b2jr + 4d^ 



kj V"lfe ) ^2j } -T W "i " -r W "j l^kj V\k , U2j ) -T W "j l^kj \ 

Ek.ibik, b2j?' + cldf Ek.ibik, b2,r + cfdf a + cf df Ek.ibik, b2jy^ + c^df 

KEkjibikM.f' + cjdf Ek.ibikM.Y' + c^df Ek.ibikM.r+cldf a + cfdf 

where 

(6ife,62,)"^^6?,6|-te, (5.6) 
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The main difference between Eq. (5.5) and Eq. (4.3) is that in Eq. (4.3) there was only one b, 
whereas here, we have a set of pairs of 6's, i.e. {(6in,&2n)}- So, for these new 6's, the first lower 
index indicates whether it comes 1st or 2nd in the cross product with v, the second lower index 
indicates what term in the sum it belongs to, and the upper index indicates which component of it 
is being taken. In our point in the middle set up, the q's consist of 5 edges and their 10 triple cross 
products; there are 15 corresponding dj's. The J2n form has a nicer appearance, but the J2ki form 
is much better when doing computations. The k subscript in jk has been omitted in the matrix 
form given. 

Using Eq. (5.4) and taking the trace and the antisymmetric part of the equation give, respec- 
tively, ^ 

4(a-l) + Y,Ci-di = (5.7) 

i 

and 



B-^C = (5.8) 



where B and C have matrix elements^ 



5"^^E^"VL-&2. (5.9) 

jk 

and 



t 



These (un-contracted) identities are in terms of deficit angles, internal triangle angles, normal- 
ized areas and normalized four-volumes. The normalized areas, naturally 2-index tensors, appear 
explicitly only in the tensor identity, Eq. (5.8). 

Using the explicit form of Ri^^tot iii Eq. (5.5) and Eq. (5.8) the completely contracted Bianchi 
identities can then be written as 

1 1 

~2W{v) f- 



J2 ea/37A (2i3 - Cr^l^ = (5.11) 



To first order in deficit angles, Eqs. (5.7) and (5.8) reduce to, respectively, = and 

B''^ = J2^^^''^uK = (5.12) 



■'In d dimensions, 4(1 — a) ^ d{l — a) 



In d dimensions, there are d — 2 &'s contracted with the Levi-Civita tensor. 
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The above equation can be rewritten, using Eq. (2.4), to give 

Y.6,uf = (5.13) 

i 

where we have substituted in for l'^ using Eq. (5.1). Following a similar derivation to that of 
Eq. (5.11), the above equation allows the fully contracted Bianchi identities to be written for small 
angles as 

-41^ E ^a^,x(T.^^Uf)l^ = (5.14) 

Eqs. (5.13), (5.14) (as well as the Eq. (8.10), reduced to four dimensions) have been verified to first 
order. For higher order, these equations have been shown to be violated in favor of, respectively, 
Eqs. (5.8) and (5.11) (as well as Eq. (8.11), again reduced to four dimensions). Also, Eq. (5.7), a 
trivial = to first order, has been shown to be true generally. ^ 

6 Discrete Riemann Tensor and its Dual 



In the previous two sections an explicit, exact form of the lattice Bianchi identity was derived 
by considering a product of rotation matrices (in three and four dimensions respectively) along 
a homotypically trivial path. Curvature enters the above lattice Bianchi identities through the 
deficit angle, but one interesting question left partially open is the relationship between the exact 
lattice Bianchi identities (which involve products of rotation matrices) and the continuum Bianchi 
identities (which involve derivatives of the Riemann tensor). In this and the following section an 
expression of the Riemann tensor and its dual will be derived, which will eventually be used to 
show that in four dimensions (Section 9), if one proceeds from the continuum Bianchi identities (in 

integrated form) and inserts the expression for the discrete Riemann tensor discussed below one 

^We note here that these results apply to the Lorentz case as well. One sets ct — ~ix4, and finds the x, t rotations 
equivalent to the x, X4 rotations. Since all the Bianchi products of rotations, when written out in terms such as the 
Ta^b gauge transformation of Section 3, end up canceling to the unit matrix via products such as Ta^b^a^B' ^^^ 
as each such Euclidean gauge transformation has an equivalent Lorentz transformation, the Regge Calculus Bianchi 
identities hold in the Lorentz case as well. In finding these equivalent transformations, it is important to note that 
any rotation plane will have at least one spatial axis. If the other (perpendicular) axis is also spatial, one has the 
standard rotation by the computed deficit angle. If the other axis is time-like, one first computes the deficit angle as 
in the Euclidean case, multiplies it by i to get the appropriate deficit angle to use in the sines and cosines, and ends up 
using hyperbolic sines and cosines. For this case, it is important not to subtract out any multiple 2tt from the original 
Euclidean deficit angle before it is multiplied by i. If the other axis is null, there is no rotation, no matter what 
the deficit angle would be in the Euclidean case, as a null vector and a spatial vector cannot rotate into each other. 
These three cases are, respectively, equivalent to A'^" A]^^, > 0, A'^" A^^, < and A'^" A]^^, — 0; alternatively, they are 

also equivalent to, respectively, Al''-''~' A';s,...s,_, < 0, ^^'^ ■'"-^A^,,...,,_3 > and Al''-''-' A';s^,„s,_, = 0. 
Lastly, it is important to note that, once one fixes a going around coordinate system, one needs Euler angles to write 
a space-time rotation. 
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obtains, after promoting infinitesimal rotations to finite rotations, the same exact lattice Bianchi 
identity discussed previously in Section 4 (in Section 10 the same set of result is presented in three 
dimension). 

Having obtained the necessary ingredients to describe arbitrary lattice rotations of vectors, 
we proceed next to derive a general form for the lattice Riemann curvature tensor. An explicit 
form for the Riemmann tensor in term of rotations will be quite useful here, since it will allow us 
to establish a relationship (known to exist in the continuum) between the Regge field equations 
and the lattice Bianchi identities derived in the previous sections in three and four dimensions. 
Furthermore it will allow us to check the overall consistency of our results, since later in the paper 
we will re-derive the lattice Bianchi identities by starting from their continuum expression in terms 
of the Riemann tensor, and by promoting infinitesimal rotations to finite rotations will be able 
to show that the resulting lattice Bianchi identities are in fact identical (in d=3 and d=4) to the 
expressions previously derived in the preceding sections. 

Consider moving a vector V once around a Voronoi loop, i.e. a loop formed by Voronoi edges 
surrounding a hinge. The change in V, denoted here by 5V, is then given by 

<5y" = (R - i)°^y'^ (6.1) 

where R is the is the rotation matrix associated with the hinge. Now, 6V in the continuum is given 
by 

5v^ = i i?v ^r y^ (6-2) 

where A^^'^ is the antisymmetric bivector representing the loop area. So we make the tentative 
identification 

ii?v^r = (R-i)"/3 (6.3) 

Now, because of the sums on ^ and i', it is not immediately clear how one can divide by Ajf 
to solve for the Riemann curvature tensor. Let us, then, simply take a frame where there are only 
two components of A^^'^ , namely A^ and A^^. In that frame, using the antisymmetry of Riemann 
in /i and v, one can divide, and one finds 
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Multiplying by unity in the form of A12 (the 12 component of the index-lowered version of A^) 
over itself, one finds 

R"^i2 = ^^^^ (6.5) 

where the area of the loop is an invariant satisfying A^ = A^^A^'^ /2. Of course, all values of the 
Riemann curvature tensor with a or /3 neither 1 or 2 are zero in this coordinate system. 
Now, the above equation can be rewritten in this coordinate system as 

(R - 1)%^L 
^ V = ^ JT^^ (6-6) 

However, this is now a tensor equation, and hence valid in all coordinate systems. Still, there 
is a basic problem with this, and that is that (R — l)o-/3 is only antisymmetric to first order (for 
example, (R — l)ii = only to first order), so that this identification of the right hand side of 
Eq. (6.6) with the Riemann curvature tensor and all of its symmetries can only be valid to first 
order. So, let us write Eq. (6.6) to first order in the deficit angle 5: 

where, in our coordinate system, which we now take to be orthonormal,® 



U, 



af3 



f -1 \ 

1 



\ y 



{61 



The general form of Ua/s, in an arbitrary coordinate system, is 

Uap = ec,^^X^ll/2Ah (6.9) 



where ^1 and I2 are the two hinge vectors and Ah = J A\^A^^ /2 with Ah the area of the hinge and 
A):^^ the associated bivector.^ Note that 5 is invariant because it is formed from angles which are 
arc cosines of ratios of invariant dot products to invariant lengths, and furthermore that eapnu = 
yjg [1234] with [1234] = +1. So Eq. (6.9) is a tensor, and hence so is Eq. (6.7). Since the expression 
in Eq. (6.7), after lowering the first index, also has the appropriate antisymmetry in the first two 
indices, we can see it as a Riemann tensor candidate. 



*A minus sign appears above in the Lorentz case of a x-t loop. 

^In the Lorentzian case of a space-time hinge, An — yJ—A\^A'^"'/2. 
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Now, by going back to our special coordinate system, it can be easily verified that the dual of 
the normalized loop area is the normalized hinge area, i.e. 



^r 2 Ah 

This allows us to rewrite Eq. (6.7) simply as 

«V = ■£ ^"pU^- (6-11) 

This expression has been discussed by the authors of Ref.[7, 10], and is known to possess all 
requisite symmetry properties of the Riemann tensor. However, as these authors note, it implies 
that the square of the Ricci scalar, the square of the Ricci tensor and the square of the Riemann 
tensor are all directly proportional, with constants of proportionality independent of the edge 
lengths. We now define 

^ V = ^ V(^' ^) ^ E ^ V(^) (6-12) 

hCb 

for a formulation of Regge Calculus which breaks up space-time into d-boxes. Here, one chooses 
one particular vertex v to be the origin for the g^u of d-box b. Then the metric can, for example, 
be defined via dot products of the box's non-diagonal edges containing v, or, equivalently, can be 
defined via edges in the box containing v which are in one of the box's d-simplices. Other vertices 
V for other d-boxes are chosen to be in the same relative location within the d-boxes, so that each 
hinge curvature term R^^^^ih) is included once. 

Two other definitions of the Riemann tensor are possible in principle, but as we will be show 
below, fail to pass a crucial test. Noting that each rotation plane has four hinges for a d-box, one 
could have alternately defined 

^V^iE^VW (6-13) 

hCs 

This could be useful in a situation with large curvatures where the increased computation due to 
averaging and, as we will see, more complex Einstein tensors GS' 's, is more than canceled out by 
being able to use fewer d-boxes due to the relaxing of the origin vertex location constraint. 

Alternatively, if one uses only d-simplices and does not require them to form d-boxes, then one 
can define 

^ V ^ ^ V(«' ^) ^ E TT ^ V(^) (6-14) 

hCs ^^"^h 

where Ng^ is the number of d-simplices containing hinge h. 

Now, one objection to these definitions is that the number of independent curvature components 
used per d-volume, for these three cases respectively (Eqs. (6.12), (6.13) and (6.14)), is d{d— l)/2, 
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2d{d— 1) and d{d+ l)/2, none of which is the expected d'^{S — 1)/12, the number of components of 
Riemann in the continuum. However, let us now calculate the number of components of Riemann 
about a vertex in the d-box case where consistent labeling of the origins are possible, i.e. where 2'^ 
d-boxes surround each vertex. The number of m-boxes surrounding a given vertex is 2"^(m) where 
we choose m of the d axis directions, and note that there are two choices per direction. Setting 
m = d — 2 gives 

Nh,v = 2'^-^d{d-l) =d^{d^ -l)/12 d = 2,3 (6.15) 

> d^{d^- 1)/12 d>4: (6.16) 

So, if one sets up a going around coordinate system at a vertex, and defines the curvature via 

R^^^s^T.R'^P'ysW (6.17) 

hDv 

the symmetries of the form of the hinge curvatures keep the number of independent components 
of the curvature equal to d'^{d'^ — 1)/12 in any dimension. This equation, then, having along with 
the natural form for the Riemann curvature tensor, all the symmetry properties of the continuum 
Riemann curvature tensor as well as the correct number of independent components, is the Regge 
analog of the Riemann curvature tensor. ^^ 

We are now in a position to calculate the dual Riemann tensor via 

j^af3^S ^ l^^uySj^aP^^ ^ ^5 e^^^'^U'"^ A^JAf (6.18) 

where the loop now used is one of the aforementioned inner-d-box loops. ^^ We can interpret the 
above result in the following way. Here [/"^j corresponds to the "effective" C/"^ of the loop, and S 
is the "effective" deficit angle of the loop. The "effective" Riemann curvature tensor is the sum 
of all the "old" Riemann curvature tensors corresponding to the areas that the 2 — 1 = 15 edges 
form, so that the effective 6 and If^^ are determined for each standard square loop (which can only 
encircle areas formed by the 15 edges). Note that, since our "old" Riemann curvature tensors had 

^''Were we to use d-simplices which need not form d-boxes, then the consistent choice of origin vertices may be 
either impossible or very difflcuh. One might also make this choice for a d-box situation with high curvatures where 
the increased computation due to averaging is more than canceled out by being able to use fewer d-boxes. But this 
is more complicated and our results for Gj'' will also show a distinct preference for the first option. Still, in cases 
where curvature is sufficiently high and computing power will not allow for a sufficient breakdown of space-time into 
d-boxes such that a consistent choice of origin vertices is possible, this definition may be useful. Indeed, for this 
reason, were one to break down space-time into d-simplices without forming d-boxes, this definition would likely be 
necessary. Still, in such a case where either a quick estimate is needed or where computing power is insufficient to 
use the first definition, this definition would then be useful. 

^^In the Lorentzian case, there is a minus sign in this definition of R"^''^ . 
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all the symmetry properties of the continuum Riemann curvature tensor, and because our "new" 
(i.e. actual) Riemann curvature tensor is the sum of these old tensors, it too has all the symmetry 
properties of the continuum Riemann curvature tensor. And, further, it has the advantage of 
allowing us to consider, in the case of small deficit angles, only perpendicular loops (more precisely, 
perpendicular loop areas), as we do in the continuum. Now, using Eq. (6.10), Eq. (6.18) can be 
written as 

which, incidentally, holds in any number of dimensions simply by adding d — 4 extra indices after 
the 5 indices. As the action will involve an integration, it will now be helpful to find out the total 
four-volume enclosed by the hinge and its loop in terms of the product of the loop area and the 
hinge area. 

In the 4-d case the denominator. Ay Ah, is six times the four- volume enclosed by the hinge and 
the path, here denoted by '^V{h,T), this latter four- volume being formed via the Voronoi vertices 
surrounding the hinge and the hinge vertices. To see this fact, note that, in the 3-d case, as one 
approaches a vertex, the area perpendicular to the hinge is similar to the loop area, but gets smaller 
proportionally to the square of the perpendicular distance from the point, labeled p, which marks 
the intersection of the loop area and the hinge. (In particular, all lengths defining the area are 
proportional to this distance, and hence areas, as functions of products of pairs of lengths, must be 
proportional to the square of this distance.) By letting lia be the distance from one of lis vertices 
to p, and letting In, be the distance from lis other vertex to p, one finds the total 3-d volume 
formed by the loop and h to be 

^V= f '\Ar){s/hafds + / '\Ar){s/hbfds = Arh/3 (6.20) 

Jo Jo 

In 4-d, the total four- volume formed by the loop and the hinge is obtained by noting that the three 

volume shrinks proportionally to the cube of the perpendicular distance from the three volume to 

the remaining vertex, so that, letting I2 be the component of the remaining hinge edge which is 

perpendicular to the tetrahedron. 



'V= / {Arhms/lifds = e^,.,sA^''llli/2A = e^.^^sA^" Af /2A = ArAh/& (6.21) 

Jo 



/o 

where the next to last term is found by taking our special coordinate system, with two axes alon^ 
Ay and two axes along Ah. In d-dimensions, one has V = £^1/75... ^p'^V^ '" /d\ = 2ArVh/d{d — 1) 
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Incidentally, the only constraint on the loop here is that it be perpendicular to the hinge. One has 
therefore 

j^afSyS = 1 J_ u^'f^Al^ (6.22) 

Using Eq. (6.11), we find 

S ea/37i-7d-2^h e '"'^7i-7;; 2 S , , 



I=- / R^d^x = -Y,l ^^d^x (6.24) 

Jx ^ Jh Ay 



so that the Einstein-Hilbert action becomes 

2(5 

h 

where X is all space and /^ is defined as the integral over '^V{h, Th), the d-volume formed by the 
hinge and the loop area. As the integrand is constant over the hinge, we do the integral by simply 
calculating this aforementioned d-volume. One can now can write the action as 

where Vh/d{d— 1) is the (d-2)-volume of the hinge, reducing to an area A^ in 4-d. (The sign in the 
above equation becomes positive in the Lorentz case.) 

7 Action Variation and Einstein Tensor 



This section will be devoted to discussing the relationship between the expression for the dual 
of the Riemann tensor, as given in the previos section, and the Regge field equations. We will show 
that the above construction is indeed consistent, by deriving from it the Regge field equations. 
While this does not constitute a general proof of correctness of the proposed expression, it does 
provide one rather significant test. 

Next one would like to vary the action with respect to g^y. In the continuum each g^i^ is an 
independent variable upon which the geometry is based, which is why one varies the action with 
respect to them. Since, in Regge Calculus, the geometry is completely based on edge lengths, only 
edge lengths can be used in computing g^y. In particular, any g^y for the loop plane is not a true 
g^jtu, since these quantities can be defined independent of edge lengths. Furthermore, one can see 
that if one did vary the above action with respect to such a quantity, one would get zero because 
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the hinge areas are independent of any quantity which is independent of the edge lengths. So, one 
is led to consider the metric (7^1^ as being defined as a function of hinge edge lengths. Here, we take 
one {d — 2) x {d — 2) gf^^, matrix per hinge. Using the invariance of 6 in any dimension[l], '^V^ as 
the Voronoi d- volume surrounding the hinge and Vh = Ji^ \/\^d '^x/{d — 2)!, we now compute, for 
an individual hinge, 

'pi' 



^n, - X^^G,^ ^\9\dx = -^(^3T)^7P^^5^^ 



la\'' 



9-yp r 1 c- rt jd—2 

[ai...ad-2\ S'<7il---5a ,_i7'-lC><xy5<7 /_^i7' + l---S'ad_2d-2 « = 



Ihdl 2^\g\ 

^ -S^^'ShVh (7.1) 



d{d-l) ^ 

To validate this result, ^^ which confirms the correctness of using (d — 2) x {d — 2) hinge metrics, 
we compute the above quantity in a different manner via^^ 

K ^^ Vl^l 2{d- 3)1 J.iv, 

^ ^^ r TTal3^rlSi...5d~3 2 

' 2{d - sy.ArVh '"''''''■■■'''-' '^ ^ " d{d-l 



p,'S....s,_.AU"^V;:''-'''' = -^^nH\^,'' ^hVu (7.2) 



Now, noting that this final result is independent of the metric, and noting that the br? indicates 
that the components of QS^ are most naturally taken along edges, we find 



h. C d — b ox 



In a given d-box of our global coordinate system, where these Ts are the box's (non-diagonal) axis 
/, the integrated Einstein tensor is defined via hinges within that d-box, so that there is no global 
over-counting of each contribution to Q^^ . We now easily see that, using our d-box's g^^^ defined by 
the axes' edges, 

^"' = -7^7:7^ E ^hyhQiv and Q^^' = -—^ ^ b^vj^' (7.4) 

d{d - 1) ,3f- did - 1) ^^f- 

/i C d — b ox /i C d — b ox 

The cosmological constant term is derived in each d-box from the d-box's g^^y via the variation of 
2 / Ay^ d X exactly as it is in the continuum to get 

\g and hence A(5; and Xgiy (7.5) 



^^The sign is opposite for a space-time loop, because Eq. (6.25) changes sign for the reason given in the footnote 
for Eq. (6.7). Also, letting ^fg — > \J~g for space-time hinges has no effect on the sign. 

^■^Once again, the sign of the resuh changes for a space-time loop for the reason given in the footnote to Eq. (6.7). 
For any Lorentzian hinge, the minus sign referenced in the footnote to Eq.(6.18) is canceled by the minus sign in the 
inverse to Eq. (6.10). 
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This is the most natural result for the left hand side of Einstein's equations for Regge Calculus 
because of its metric independence. A computation of Q^'^ or t/^jy for edges would have been 
counter-intuitive because one would be using a different metric for each hinge in the sum. Also, 
using Gf^S' would have explicitly involved the areas of the hinges' loops, which are not uniquely 
defined. Now, the standard result for the field equations in four dimensions 

- ^ (5fe cot Oop (7.6) 

is derived by varying the integrated action with respect to the edge lengths, and hence is most 
similar to Q^^^ , which is found by varying with respect to g^^, which is a function of squares of 
edge lengths. Hence this standard result is not metric independent. However, as long as one has 
rp^iv _ Q^ there is no problem with it as Regge Calculus then only involves edge lengths. However, 
if one wishes to include T^^ , and if this T^^ is not derivable from an action, then one must use 
Eq. (7.3). Also, Eq. (7.3) seems more natural, and, except for having a more limited sum, has the 
same form as the action. Lastly, we note that the individual terms in Eqs.(7.3) and (7.6) can be 
related via the transformations gn = if, (722 = ^2 ^^*^ 5'12 = (^1 + ^2 ~ ^^)/2 so that the variation 
of the integral of the action can be transformed, for any one fixed term, between the two sets of 
variables via a Jacobian. 

Ideally, one wants a one-to-one mapping between the (independent) edge length variables and 
the independent g^^i, variables. This is easily achieved by choosing as the independent ^^j^'s the 
g^fi's, i.e. the metric components equal to the square of a given edge length. Of course, we can 
always switch global metric variables by replacing some ^^^'s with an equal number of g^u's (while 
being sure to maintain the relative independence of the new set of metric components). We do this 
now, choosing gn = if for all [d — 2) x {d — 2) metrics for the hinges bordering /i; we let all such 
metrics be based at one of /I's vertices. All other lattice g^u^s are then chosen such that the set of 
all g^u^s are independent. A simple way to make such a choice is to choose the metric components 
equal to the edge length squared of each edge not in a simplex bordering li. Now, replacing /i with 
L and varying with respect to gn = L? gives, for 4-d 

which is easily related to the standard result by multiplying by —QL. 

As an exercise we show next the equivalence of our Einstein tensor to the variation of the action 
with respect to the edge lengths in the case of four dimensions. For a hinge with edges li and I2, 
taking gu = if, g22 = I2 and gi2 = h-h = ('1 + '2 ~ ^^)/2 where L is the third edge of the triangle 
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formed by /i and I2, we have 



A = -\/giig22 -912921 (7.8) 



where, for variation purposes, we take gi2 and 5(21 to be independent. This imphes 

dA _ 522 dA _ gu 

dgii ~ 8.4 9<722 ~ 8^ 

dA g2i dA 512 



(7.9) 
(7.10) 



dgi2 8A dg2i 8A 

and therefore 

dA dA dg^u gi2 L 

where we have used gi2 = hh cosOop, A = I1I2 sin9op/2, gi2 = (/f + ^2 — -L^)/2, gu = ll and (722 = ^i- 
This is consistent with the standard Regge Calculus result for the field side of the Einstein field 
equations. Also, similar results apply for either of the other two edges if we change the metric so 
that its diagonal components do not have that edge. If we do not change the metric, we find, for 

e.g., for /i, 

dA 1 

—- = -(^2 CSC 6'op- /i cot 6iop) (7.12) 

OL\ 4 

where 9op is still the angle opposite edge L. Also, one can "go the other way", and calculate ^— 

via 

dA dAdLdL^ (L ,. Wl^n cot^op 

which is consistent with Eq. (7.10), where we have used L^ = gu +522 —512 —521 as well as /f = gu 
and I2 = 522- So, for each edge corresponding to a G''' ,, it is easiest to vary each of the attached 
hinges with respect to the metric none of whose diagonal elements corresponds to that edge. 

Next we show that G^ , is proportional to A 5 for 7 = 7', and zero otherwise, independently 
of which edges were chosen to form the diagonal elements of the metric. So, the contribution from 
each hinge to the Einstein tensor for an edge is independent of the form of the metric if one is 
considering the one index up and one index down form of the Einstein tensor. So, Regge Calculus, to 
be "metric independent" , would choose this form of the Einstein tensor to work with. In particular, 
this metric independence does not hold when considering G^jy or G^^^ . With a metric independent 
form one can have some sense of justification when one adds the various hinge contributions to edge 
L's Gj^; G and Gll, not having metric independence, are undefined, or, at least, non-unique. 
Interestingly, the variation of L contributes only to G'^'^ and G'^'^, but not G^^] G^^ receives 
contributions from varying li and 12- In terms of the number of G'^^'s, we are consistent: there 
would be three gpaS per triangle, giving three G'^^'s per triangle, which we take to be G ^^ , G \ 
and G^^. 
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A crucial question at this point is whether the form of the Einstein tensor, Lcot 9op/2, is metric 
independent. This form is obtained by varying with respect to edge length L, which is similar to 
varying with respect to edge length L^ or li-l2. Hence, this form most closely corresponds to either 
G or G^^, and therefore is not metric independent. This can also be seen by noting our previous 
result of Eq. (7.12). 



Next we calculate G \ and obtain 
'1 



11 I „ ^12 



G\ =giiG''+gi2G 



2 I 922\ , f. J a \ f-hh COS 9op 



/?(||)+ai/2COsM( 



8^ 



2(5 = A6 (7.14) 



which shows that our procedure is consistent. 

8 General Form of the Bianchi Identity 

In the previous sections on three-dimensional (Section 4) and four-dimensional (Section 5) 
rotations an exact and explicit form of the lattice Bianchi identity was given in terms of product of 
rotations along homotopically trivial paths. Later on in the paper an expression for the Riemann 
tensor and its dual were proposed, which among other properties, correctly reproduce the Regge 
field equations. In this section we will show that if one proceeds from the continuum Bianchi 
identities and uses the above given expression for the lattice Riemann tensor then the following 
is true. Firstly, the resulting lattice Bianchi identities coincide with their weak field counterparts 
discussed in the sections on 3-d and 4-d rotations. Secondly, when the infinitesimal form of the 
rotation matrix is promoted to the correct finite rotation expression (of Sections 4 and 5), the 
same exact identities derived previously are obtained. These result presented here are therefore 
intended to bring out one more time the close relationship between the Riemann tensor (in terms of 
which the continuum identities are naturally formulated) and the finite rotation matrices (in terms 
of which the exact lattice identities are formulated). A lattice expression for the Riemann tensor 
provides a natural bridge between these two different realms. 

We will start here by considering the lattice Einstein tensor. Since the Einstein tensor's com- 
ponents live on edges, derivatives of the Einstein tensor must occur at vertices. Since the Einstein 
tensor is discontinuous at vertices, the easiest procedure is to integrate the Bianchi identities over 
a d- volume around a vertex and divide by that d- volume, at least in the case of small deficit 

angles. ^^ Then, the result is generalized to arbitrary deficit angles. So, we begin by deriving the 

^■'The assumption of a flat coordinate system, so that edges have definite directions and so that we wiU be able to 
use Gauss' theorem, introduces errors of 0{5 ). 
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un-contracted, the partially contracted, and the completely contracted Bianchi identities in terms 
of the dual Riemann tensor, which will be the dual of the local Riemann curvature tensor as op- 
posed to the global form, for reasons discussed earlier. In the continuum the un-contracted Bianchi 
identities read 



-"- a/3;7 "•" -"" /37;a "•" ""■ 70;/? ~ '-' 



(8.1) 



The partially contracted and completely contracted Bianchi identities are obtained via contraction 
of Q with a' or /and /? with /3'. Using the d-dimensional dual Riemann tensor, one can write the 
un-contracted Bianchi identities into a divergence form via 






■ i'd-2''^ );a^\'^^aai...ad-2-'^ /;/3J ~ '-' 

(8.2) 



(8.3) 
^^^(e„/375....5,_3^"'^'''"-'^-^);A = (8.4) 

Let us first consider the case of small deficit angles. The next step is (1) to integrate over and 
divide by '^V{v), the d- volume surrounding vertex v, (2) to use Gauss' theorem, and (3) break 
up the (d-l)-dimensional surface integral into (d-l)-dimensional surfaces composed of (a) the 2-d 
Voronoi rotation areas and (b) the (d-3)-dimensional portions of the hinges' surface volumes (the 
'^-^Vh's). One finds 

1 



{d-sy.'^vii 



E/ 



al3jSi...Sd-:i 



R 



:a'l3'X5i...5d-3 



fix \/\g\ d'^-^x = 



(8.5) 



Rewriting R using Eq. (6.7) and a generalized form of Eq. (6.10), letting the index A -^ Aq, and 
defining the volume tensor 

1 



TrpOPl---Pd-3 



id -2)1 



E 



V ^) ''h,io''h,h---'-h,id- 



-3 
id-'i 



(8.6) 



(io,Ji,...,Jd-3)=7'{0,l,...,d-3) 



where (—1) is the sign of the permutation, we find the integrand to be 



-Q!/375i...<5d_3 



Q'/3'd-2-(/-'*o5i...5d_3 
h 



5h U'^'f''-'^V, 



Ar '^-^Vh 



e^iyAoAi...Ad-3^r^(^/i,l — ^/i,o) ^■■■i^h,d-3 — Ihfi) "^ ^ 



(8.7) 



Note that for each hinge the above expression is independent of the choice of assignment of the 
labels //i^O) ^/i,i; •••) ^h,d-z to particular edges of the hinge. Furthermore when Eq. (8.6) is substituted 
into the above equation, each term in Eq. (8.6)'s sum can be replaced with +Iq^Ii ■■dj'Z^ , so that 
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the contraction of edges with the 2-nd Levi-Civita tensor always gives iAr Vh where (4) we 
define A^ such that the minus sign is always chosen (so as to be consistent with counter-clockwise 
rotations in 3-d). After canceling the integration with AY'^~^Vh, Eq. (8.5) can then be rewritten as 



ea/375i...5d_3 



E'^/^^"'''Efc-& = (8.8) 



where hb refers to a "hinge base", and is a (d-3)-simplex, and hence necessarily a "base" to the 
hinges containing it. One can then invert the sums to find 

- (^^t(v) ^^^'^ <) ^ti-te-3 = (8.9) 

where a unique labeling assignment has been chosen for each hinge base, Here U^l^ = itt/" " , the 
plus or minus sign chosen when {lhb,i7 ■■■i^hi,,d-3} is, respectively, an even or odd permutation of 
{lh,i,-,lh,d~3}- I-e-, the sign of U^f^ is chosen such that ea/37<Si...5d_3^^,fe,^o^^j„i---V67i-3 < ^r 
one, and hence for all via projection to 3-d, h D hi,. ^^ 

Using as the equation for V^^^ " ''"^ a (d-3)-dimensional version of Eq. (8.6), the above equation 
can then be rewritten as 



^ ^ htDv hDhi, 

which for the special cases of three and four dimensions reduces to the approximate small deficit 
angle expressions discussed previously in the sections on 3d (Section 4) and 4d (Section 5) rotations. 
Using J2^h U'^'^' ~ E{Rot - 1)"''^' « {n[l + {Rot - 1)] - 1}"'^' = [U{Rot) - If'^', where "f«" 
symbolizes equalities to first order in deficit angles, the above equation generalizes to, for arbitrary 
deficit angles, ^^ 

- '^X(T' ^ [ n (Rot,,,^) - 1]"'"' <■■"'-' = (8-11) 

where the product of rotations is over a null path, and where /?' and a' refer to, respectively, the 
first index of the first rotation matrix and the last index of the last rotation matrix. In the special 
case of three and four dimensions the above result is equivalent to the exact expressions discussed 
previously in the sections on 3d (Section 4) and 4d (Section 5) rotations. ""^^ 

"'^Here, we replace A"^ '^ with A'^ j^ (x +U^ /^ . 

^^Here, ]^(-Rot) is a product of mixed tenors with the first index up and the second down; however the first index 
of the last rotation matrix (on the left), a' , and the 2nd index of the first rotation matrix, /?', can be either up or 
down. 

^^In the Lorentzian case, the above equations have a plus sign; each term of the above equation gets an additional 
minus sign either because of the reason noted in the footnote to Eq. (6.7) or because the presence of the time 
component in a hinge would otherwise cause a reversal of the direction (and sign) of A" f^ due to the change in sign 
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Let us add a few comments regarding the result just obtained. The coordinate system in which 
this equation is written is one of the possible global coordinate systems described in the section on 
Geometric Setups (Section 3). By taking the Levi-Civita tensor out of the J2htDv^ '^^ ^^^ taking 
the Levi-Civita tensor to be defined at vertices such as v}^ Now, one can ask how we can use a 
d-dimensional Levi-Civita tensor when in our derivation of ^^ we only used {d — 2) x {d — 2) metrics; 
the answer is that there we were considering Voronoi volumes which had two dimensions not set 
by edge lengths, whereas here all dimensions of every volume emanating from vertex v involve edge 
lengths. 

Also, the partially and fully contracted Bianchi identities, now for arbitrary deficit angles, are 
easily obtained via contractions. The above Bianchi identities can be converted to scalar form 
via dot products with d d-vectors which form a d-volume in our going around coordinate system 
discussed in Section 3. Do note that, when one changes going around coordinate systems, dot 
products do change, so that even here, what is a Bianchi identity for one going around coordinate 
system may not be for another going around coordinate system. Another (going around) coordinate 
system difference would be that F's such as that in Eq. (3.4) would be unity for some coordinate 
systems, and not for others, so that a different order of rotations, utilizing the r(s) which are 
unity, would be appropriate for these other coordinate systems. Still, once the dot products are 
taken, the Bianchi identities are in terms of edge lengths, which are independent of coordinate 
system; hence, once dot products are taken, any Bianchi identity is true independent of the going 
around coordinate system. In summary, each going around coordinate system has its own Bianchi 
identities, which themselves can be projected into coordinate system independent equalities. Of 
course, once one substitutes in for the deficit angles their explicit expression in terms of edge 
lengths, one attains zero identically. However, if one does not do this, these identities give relations 
between the deficit angles and squared edge lengths. 

9 Bianchi Identities in Four Dimensions 



The discussion of the previous section was for a general dimension d. In this section we will 
focus on the four-dimensional case, and proceed from the continuum Bianchi identities to derive the 
weak field expression for the lattice Bianchi identities. Once the resulting infinitesimal rotations 
are promoted to finite lattice rotations, the resulting identities coincide with the exact form given 



18 



We needn't have done this; however, taking it out of X^^-,^ or J^^-,^ avoids unnecessary complications 
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previously (Section 5). At the end of this section we wih then coUect the exphcit formulae for the 
exact un-contracted, partially contracted, and fully contracted lattice Bianchi identity, valid for 
arbitrary manifolds. 

In the four-dimensional case the un-contracted identity in the continuum reads 

\[{eap^.uR'''^'n■a + (^Pl^^^R'''^'n,c. + {^^a^.uR'''^'n■A = (9-1) 

Now, defining the index 5 such that [a/37(5] 7^ (for d > 4, we choose one set of indices, an ordered 
set, 5 = 6i62---Sii-3, which satisfies [af3^6] / 0), we find 

{eafS-ySR'''^'^')n + (e/37a5i?"'^'"');a + (6^^/35 i?"'^'^');/3 = (9.2) 

where, as we have stated, a, /3, 7 and 6 are not summed anywhere in the equation. We can now 
rewrite the above equation as 

{eat3y5R'''^'^');X = (9.3) 

where A is summed over, and 6 is not. Actually, in 4-d, the above equation does not change even if 
6 is summed over all directions; for d > 4, one simply multiplies by a factor of l/{d — 3)! to get the 
exact correspondence to the un-contracted Bianchi identities if one decides to sum on 5 (which, for 
d > 4, becomes 5, and the "sum on 6" is the sum over all the J's of 6, which has non-zero terms 
only when none of the J's are equal to a, /? or 7). 

To get the partially contracted Bianchi identites, one simply either sets a = a' or /3 = /3', and 
then sums over either a or (3. Setting (3 = jS' gives 

(ea/375^"'''^');A = (9.4) 

Here, we use the Einstein summation convention for all indices, so that for d > 4, we would need 
a factor of l/(d — 3)! from Eq. (9.3). In particular, we note that 5 must be summed over, because 
we are summing over two values of /3, and therefore 5 must run over different indices (or, for d > 4, 
different /3's implies that not all (5's which give non-zero terms will be comprised of the same set 
of indices). To get the fully contracted Bianchi identities, we set a = a' and sum over a in the 
partially contracted Bianchi identities to find 

(e,^^,i?"^^'^);A = (9.5) 

where, once again, we would have the factor of l/(d — 3)! for d > 4. 

The next step is to interpret the un-contracted Bianchi identities, Eq. (9.3) in Regge Calculus. 
Using the lattice expression for the dual of the Riemann curvature tensor, Eq. (6.22), we have 
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:(^4 U"'^'A^%, = (9.6) 



We now integrate this over a vertex v, and divide by the sum of ah Voronoi four- volumes surrounding 
that vertex, ^I^. 

One can rewrite this in terms of edges, as the components of the Einstein tensor hve on edges, 
and because the completely contracted Bianchi identities are written in terms of components of the 
Einstein tensor. 

^ 4^ E / (^ ^^ ^"'^' Af );A ^\d'x = (9.8) 

Since each hinge appears on two edges, we find 

^if EE/(^ 5^ U-'''An,,^\d^. = (9.9) 

which implies that the dual Riemann tensor for an edge is one half of the sum of the dual Riemann 
tensor for all hinges containing that edge. Using Gauss' theorem, one then finds 

^I^EE /(^ ^^ U-'^'^n) n,^\d'x = (9.10) 

Noting that the radial direction is contained within the hinge, one can integrate over the loop area, 
leaving the only integration variable to be the direction ib(L — /), where L is the other edge of the 
hinge which also contains (and is directed outward from) the vertex v. L — I here is the unique 
direction, up to a sign, perpendicular to the surface normal without any component contained 
within the loop area. We choose +(L — I) to ensure consistency throughout the sums. One has 
therefore 

^i^EE /(^r^ 4 U-'P'Al') n^^\ds = (9.11) 

Noting that A'^ = L^l — l^L (where we have chosen a sign convention for A^), one can write the 
dot product L ■ h as 



" Ar\L-l\ ~ Ar\L-l\ ~ \L - 1\ ^^'^^^ 



Then using this result and a similar one for I ■ n one obtains 
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\wT.i:'-^5>^U-'^'ArAUl-Lf = (9.13) 

which is equivalent to 

^ Y, e^p,,{±) Y: 6, U^'P'il -Lf = (9.14) 

" iDv hDl 

where the it sign will indicate whether the hinges, for each particular hinge base, are gone around 
in a counter-clockwise or clockwise fashion. This is then further simplified when one considers that 
the contribution for each edge is the same whether it be in its L (non-hinge base) contribution, or 
its / (hinge base) contribution, noting that I ^^ L leads to it ^ =F. Indeed, given that we started 
with a sum over hinges, Eq. (9.8), and broke that down into a sum over edges, Eq. (9.9), each edge's 
contribution to the hinge must be identical. So, we can take each edge's / contribution, double it, 
and omit its L contribution to find for small deficit angle 

^ Y e^p,s{±) E ^h U^'^'l' = (9.15) 

^ iDv hDl 

The partially contracted and completely contracted forms are then easy to obtain. The partially 
contracted form can be written for small deficit angle as either 

iTT E E e"/57^ '^^ U^'^il^ -I)' = (9.16) 



'" iDv hDl 



or 

1 



EE^"/37<5 Sh U'^^'iL -I)' = (9.17) 



4T/ 

'" iDv hDl 

and the fully contracted form for small deficit angle is 

4r E E 'c.i3ys Sh U^^{L -If = (9.18) 

^ iDvhDl 

We note that each fixed I term in ^^ can be taken to correspond to one term in (what was) the 
divergence, because each edge length corresponds to one component of the Einstein tensor G\,. 

Next we note that these small deficit angle forms can all be rewritten in terms of rotation 
matrices. For example, the un-contracted form can be written as 

^Y.Y.'-pAR<^th-ir'^'{L-if = (9.19) 

^ iDvhDl 

and for arbitrary deficit angles, this generalizes to 

^Y.'c.pAX{iRotn)-ir'^\L-lf = (9.20) 



'" iDv hDl 
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where f3' refers to the second index of the first rotation matrix, and a' refers to the first index of 
the last rotation matrix. The above expression coincides with the exact result obtained previously 
in the section on finite rotations in four dimensions (Section 5). The partially contracted Bianchi 
identities can then be written for arbitrary deficit angles as either 

j^T.'-(S,s[U(RotH) - ir'^{L -If = (9.21) 

or 

^E'-P^siUiRot,) - ir^'iL -If = (9.22) 

'" iDv hDl 

depending on the choice of indices. Finally the completely contracted Bianchi identities can be 
written for arbitrary deficit angle as 

4rE^"/375[n(^°*^)-l]"''(^-0' = (9.23) 

We can see that this has the correct units, 1/{^V), by taking the metric, and hence the Levi=Civita 
tensor, to be dimensionless. 



10 Bianchi Identities in Three Dimensions 



Let us repeat the above construction in three dimensions, without going through the details 
which very much parallel what was done in the four dimensional case just discussed. For d = 3, 
since the hinges are the edges, we get for the completely un-contracted Bianchi identities, valid for 
arbitrary deficit angles, 

3^ ea^, [lliRoti) - 1]"'''' = (10.1) 

^'' I 

Then the 3-d partially contracted Bianchi identites for arbitrary deficit angles can be written as 

either 

3^ eap, lUiR^ti) - if^ = (10.2) 

^^ I 

or 



3^ e„^^ illiRoti) - 1]°^' = (10.3) 

and finally the 3-d completely contracted Bianchi identities for arbitrary deficit angles read 

3^ eafs, lUiRoti) - ir" = (10.4) 



^V, 



38 



These expressions coincide with the exact result obtained previously in the section on finite rotations 
in three dimensions (Section 4). In three dimensions and for small deficit angles these results reduce 
to 

^Y.^aM^iU'^'''' = (10.5) 

The 3-d partially contracted Bianchi identites can be then be written, again for small deficit angles, 
as either 



1 



^Y.'c.M^iU'''^ = (10.6) 



or 



^Y.'c.M^iU''^' = (10-7) 

depending on the choice of indices. Finally for small deficit angles the 3-d completely contracted 
Bianchi identities are 

^^ IDV 



11 Bianchi Identities in Dimensions Greater than Four 



For dimensions d > 4, our 4 — d results are generalized after recalling that, in the original 
integration, each hinge appears once, so that each hinge should appear once in our results as well. 
For each hinge, all of its (d — 2) edges bordering vertex v appear in J2i once, so we need to divide 
by {d — 2). Also, for a fixed edge contained within a hinge, there are , \ = d — 3 hinge 

bases ((d-3)-dimensional simplexes which are labeled by hh), which are summed over which that 
hinge contains. So, we must also divide by (d — 3). For convenience, the order of Li,L2, ...^L^^^ 
can be chosen such that e „ pRot'^ P (L — /) > 0. Therefore for d > 4 the un-contracted Bianchi 
identities for small deficit angles read 

^ ^^ ' " iDv hbDl h=(hb,L) 

(11.1) 

where the last sum can be viewed either as over h D h^ or over edges L which, along with hf,, form 
a hinge. Taking the latter view, we can note that, for a fixed term in this triple sum, we should 
get the same term when / ^^ L. This, indeed, is the reason we were able to divide by (d — 2) in 
Eq. (11.1). So one can simply double the above result and let (l — L)'^-'^ -^ ^d~3^ which dramatically 
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simplifies the answer to 

(d 2)(d 3) W^T.^a,,si E 6,U-^tlp...lt-:i'^-'^ = (11.2) 

[a A)[a 6) Vvi-^^^^-^i h=(h„L) 

Now, as opposed to -rjzj^rjz^Yl,iDvYl,hf,Di '^^ could have written just -^:j2^hi,^v^ but this would 
not explicitly involve the edge lengths. And, a sum over edge lengths is most directly related to 
a divergence because, as mentioned earlier, each fixed / term in Yli can be taken to correspond to 
one term in (what was) the divergence, because each edge length corresponds to one component of 
the Einstein tensor. Additionally, the identities would not look as neat, as there would not be a 
natural edge to subtract from the other hinge edges to get a hinge's {d — 3)-surface volume. 

The above formulas were for small deficit angles. For arbitrary deficit angles, we find the 
un-contracted Bianchi identities to be 

(^^^)Wi)^S,g/-^l<n «-.)-il-U-')i. ^ ,n.3, 

The partially contracted Bianchi identities can be written as either 

^ ^ T.T. 'aByS [( n R^th) - if^iL - l)t, = (11.4) 



(d-2)(d-3)'^V;f-/-- "/375 

^ '^ ' " l^v hi,Dl hDhi, 



or 

1 1 



depending on the choice of indices, and finally the completely contracted Bianchi identities can be 
written as 

^ ^ E E 'aBys [( n R^^h) - ir^iL - oL. = (11-6) 



{d-2){d-3)^V.^^^j^^-^^^ .,., 



12 Conclusions 



In this paper we have derived an exact form for the Bianchi identities in simplicial gravity. In 
four dimensions these are given by Eqs. (9.20), (9.21) and (9.23) for the un-contracted, partially 
contracted and fully contracted form respectively. In three dimensions the corresponding expres- 
sions are given by Eqs. (10.1), (10.2) and (10.4), while above four dimension the corresponding 
general results are in Eqs. (11.3), (11.4) and (11.6), with an alternative but equivalent form of 
the un-contracted Bianchi identity given in Eq. (8.11). The explicit form of the rotation matrices 
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appearing in the above-quoted exact Bianchi identities was constructed explicitly and presented in 
Sections 4 (three dimensional case) and 5 (four dimensional case). 

While fairly unwieldy in their explicit form, these identities can be shown to reduce to their 
known weak field expression in the limit of small curvatures. They provide an explicit, local 
relationship between deficit angles belonging to neighboring simplices. Their existence can be 
viewed as a consequence of the local invariance of the Regge action under small gauge deformations 
of edge lengths emanating from a vertex, just as the continuum Bianchi identity can be derived 
from the local gauge invariance of the gravitational action. 

The relationship between the lattice Bianchi identities and the Regge lattice equations of mo- 
tion has been investigated as well. In the continuum the contracted Bianchi identities ensure the 
consistency of the gravitational field equations. One would expect that the same should be true 
on the lattice, in the sense that a lattice "covariant divergence" of the lattice field equations would 
identically yield zero, as a consequence of the lattice Bianchi identity. We have shown that that is 
indeed the case in the lattice theory. 

When in the Regge lattice case we go from the discrete Einstein-Hilbert action and its equations 
of motion to the fully contracted Bianchi identities, and then write down an explicit form for the 
partially contracted and un-contracted Bianchi identities as well, we find that the relation between 
the three types of Bianchi identities in the Regge theory is basically simply the contraction of 
indices, as in the continuum. Also, as the product of rotation matrices around a null path is 
critical in understanding the form of all three types of Bianchi identities, we have provided both 
a quick way to calculate this product and an understanding of, in particular, its antisymmetric 
components. Finally by appropriately projecting the lattice Bianchi identities, we have derived 
explicit expressions depending on edge lengths squared only. 

All of this is exact, valid for arbitrary deficit angles, though we have provided first order 
approximations to our results for the product of rotation matrices. An exact form for the lattice 
Bianchi identity should be useful in a variety of context, including numerical schemes for classical 
and quantum gravity. In the classical case, the accuracy of four-dimensional time evolution codes 
could be checked by evaluating the Bianchi identity along a time evolved trajectory. 



Acknowledgements 

The authors wish to thank J. Hartle and R.M. Williams for very useful discussions on the 
subject of this paper. 

41 



References 

[1] T. Regge, Nuovo Cimento 19 558 (1961). 

[2] Wheeler, J. A., 'Geonietrodynaniics and the Issue of the Final State.'. In Relativity, Groups, 
and Topology (Ed. C. DeWitt and B. S. DeWitt). New York: Gordon and Breach, pp. 467-500, 
1964. 

[3] M. Rocek and R. M. Williams, Phys. Lett. 104B 31 (1981); Z. Phys. C21 371 (1984). 

[4] M. Rocek and R. M. Williams, 'Introduction to quantum Regge calculus', in Quantum Struc- 
ture of Space and Time, edited by M. J. Duff and C. J. Isham (Cambridge University Press, 
1982). 

[5] J. Cheeger, W. Miiller and R. Schrader, Comm. Math. Phys. 92 405 (1984). 

[6] J. Cheeger, W. Miiller and R. Schrader, in 'Unified Theories Of Elementary Particles', Heisen- 
berg Symposium, Miinchen 1981, Springer Lecture Notes in Physics, P. Breitlohner and H.P. 
Diirr eds.. New York 1982. 

[7] H. W. Hamber and R. M. Williams, Nucl. Phys. B248 392 (1984); B260 747 (1985); Phys. 
Lett. 157B 368 (1985); Nucl. Phys. B267 482 (1986); B269 712 (1986). 

[8] J. B. Hartle, J. Math. Phys. 26 (1985) 804; 27 (1985) 287; 30 (1989) 452. 

[9] J. B. Hartle, in the Proceedings of the Third Moscow Seminar on Quantum Gravity, M. 
Markov, V. Berezin and V. Frolov eds.. World Scientific 1985. 

[10] H. W. Hamber, 'Simplicial Quantum Gravity', in Critical Phenomena, Random Systems, 
Gauge Theories, 1984 Les Houches Summer School, Session XLIH, edited by K. Osterwalder 
and R. Stora (North Holland, Amsterdam, 1986); Nucl. Phys. Proc. Suppl. 25A (1992) 150- 
175. 

[11] R. Friedberg and T. D. Lee, Nucl. Phys. B242 (1984) 145; 

G. Feinberg, R. Friedberg, T. D. Lee and H. C. Ren, Nucl. Phys. B245 (1984) 343. 

[12] T. D. Lee, in 'Discrete Mechanics', 1983 Erice International School of Subnuclear Physics, vol. 
21 (Plenum Press, New York 1985). 

[13] H. W. Hamber, in 'Probabilistic Methods in Quantum Field Theory and Quantum Gravity', 
(Cargese 1989, published by Plenum Press, New York 1990) pp. 243-257. 

42 



[14] R. M. Williams, 'Discrete Quantum Gravity: the Regge Calculus Approach', Int. J. Mod. 
Phys. B 6, 2097-2108, 1992. 

[15] J. Frohlich, 'Regge Calculus and Discretized Gravitational Functional Integrals', I.H.E.S. 
preprint 1981, published in Non-Perturbative Quantum Field Theory: Mathematical Aspects 
and Applications, Selected Papers. Singapore: World Scientific, pp. 523-545, 1992. 

[16] T. Regge and R. M. Williams, J. Math. Phys. 41, 3964 (2000) [gr-qc/0012035]. 

[17] R. M. Williams, Nucl. Phys. Proc. Suppl. 57, 73 (1997) [gr-qc/9702006]. 

[18] R. P. Feynman, Lectures on Gravitation, edited by F. Morinigo and W. Wagner, (Caltech 
lecture notes, 1962). 

[19] S. Weinberg, 'Gravitation and Cosmology', (Wiley, New York 1972). 

[20] C. W. Misner, K. S. Thorne and J. A. Wheeler, 'Gravitation', Freeman, San Francisco 1973. 

[21] R. M. Wald, 'General Relativity', University of Chicago Press, 1984. 

[22] W. Miller, Found, of Phys., Vol 16, No. 2, (1986); Class. Quant. Grav. 14 (1997) L199-L204. 

[23] P. Morse, Journal Class. Quant. Grav. 9 (1992) 2489. 

[24] C. Y. Wong, J. Math. Phys. 12 (1971) 70; 

P. A. Collins and R. M. Williams, Phys. Rev. D5 (1972) 1908; D7 (1973) 965; DIO (1974) 

353; 

R. M. Williams and G. F. R. Ellis, Gen. Rel. Grav. 13 (1981) 361. 

[25] R. Sorkin, Phys. Rev. D 12, 385-396, 1975; J. Math. Phys. 16 (1975) 2432 and 19 (1978) 
1800. 

[26] A. Gentle, D. Holz, W. A. Miller, and J. A. Wheeler, Class. Quant. Grav. 16 (1999) 1979; 
S. Chakrabarti, A. Gentle, A. Kheyfets and W. A. Miller, Class. Quant. Grav. 16 (1999) 2381; 
W. A. Miller, Class. Quant. Grav. 14 (1997) L199; 
A. Gentle and W. A. Miller, Class. Quant. Grav. 15 (1998) 389; 

J. W. Barrett, M. Galassi, W. A. Miller, R. D. Sorkin, P. A. Tuckey and R. M. Wilhams, Int. 
J. Theor. Phys. 36 (1997) 815. 

[27] R.M. Wilhams and P.A. Tuckey, Class. Quant. Grav. 9 (1992) 1409; 
R.M. Wilhams, Int. J. Mod. Phys. B6 (1992) 2097. 

43 



[28] M. Rocek and R. M. Williams, Class. Quant. Grav. 2, 701 (1985). 

[29] B. Berg, Phys. Rev. Lett. 55 904 (1985); Phys. Lett. B176 39 (1986). 

[30] R. M. Williams, Class. Quant. Crav. 3, 853 (1986). 

[31] T. Piran and R. M. Williams, Phys. Rev. D 33, 1622 (1986). 

[32] A. Jevicki and M. Ninomiya, Phys. Lett. B 150, 115 (1985); Phys. Rev. D 33, 1634 (1986). 

[33] M. Bander, Phys. Rev. Lett. 57, 1825 (1986). 

[34] M. Rocek and R. M. Williams, Phys. Lett. B273, 95 (1991). 

[35] H.W. Hamber and R.M. Williams, Phys. Rev. D47 (1993) 510. 

[36] G. Immirzi, Class. Quant. Grav. 11, 1971 (1994) [gr-qc/9402004] ; Class. Quant. Grav. 13, 
2385 (1996) [gr-qc/9512040]; Nucl. Phys. Proc. Suppl. 57, 65 (1997) [gr-qc/9701052]. 

[37] J. B. Hartle, W. A. Miller and R. M. Williams, Class. Quant. Grav. 14, 2137 (1997) [gr- 
qc/9609028]. 

[38] H. W. Hamber and R. M. Williams, Nucl. Phys. B487 345 (1997); B451 305 (1995); Phys. 
Rev. D59 064014 (1999). 

[39] H. Coxeter, 'Regular Polytopes', Methuen and Co. Ltd., London, 1948. 

[40] H. Coxeter, 'Regular Complex Polytopes', Cambridge University press, 1974. 

[41] P.S. Alexandroff, 'Combinatorial Topology', Graylock Press, Rochester 1956. 

[42] A.D. Alexandroff, 'Konvexe Polyeder', Akademie Verlag, Berlin 1958. 

[43] W. Thurston, The Geometry and Topolgy of Three-Manifolds, Princeton University Press 
(1978); MSRI March 2002 (http://www.msri.org/publications/books/gt3m/PDF/index.pdf). 

[44] L. Kirby, Problems in Low -Dimensional Topology, UC Berkley preprint Dec. 1995 (unpub- 
lished) . 

[45] N. Warner, Proc. Roy. Soc. A383 (1982) 359. 

[46] J. Cheeger, J. Diff. Geom. 18 (1983) 575. 



44 



[47] L. Castellani, R. Catenacci, M. Debernardi and C. Pagani, Int. J. Mod. Phys. A 19, 1961 
(2004) [math-ph/02 11008]. 



[48; 
[49; 

[50 
[51 

[52 
[53; 

[54 
[55 



F. Archer and R. M. Williams, Phys. Lett. B 273, 438 (1991). 

J. W. Barrett, M. Rocek and R. M. Williams, Class. Quant. Grav. 16, 1373 (1999) [gr- 
qc/9710056]. 

J. W. Barrett and R. M. Williams, Class. Quant. Crav. 5, 1543 (1988). 

G. Arcioni, M. Carfora, A. Marzuoli and M. O'Loughlin, Nud. Phys. B 619, 690 (2001) 
[hep-th/0107112]. 

M. Carfora, C. Dappiaggi and A. Marzuoli, gr-qc/0310039. 

V. M. Khatsymovsky, Phys. Lett. B 586, 411 (2004); ibid. B 567, 288 (2003); B 560, 245 
(2003); B 530, 251 (2002). 

G. Batrouni, Nucl. Phys. B208 (1982) 467. 

G. Batrouni and M. B. Halpern, Phys. Rev. D30 (1984) 1782. 



45 



